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AN EXPERIMENTAL INVESTIGATION OF THE FLUTTER OF SEVERAL 
WINGS OF VARYING ASPECT RATIO, DENSITY, AND THICKNESS 
RATIO AT MACH NUMBERS FROM 0.60 TO 1.10 


By Raymond Herrera and Robert H. Barnes 
SUMMARY 


An experimental investigation has been conducted to determine the 
flutter characteristics of several rectangular wings of variable aspect 
ratio at Mach numbers from 0.60 to 1.10. The wings were solid aluminum 
or steel and were 2, 4, and 6 percent thick. 


Flutter was encountered at both low and high angles of attack. At 
low angles of attack, the flutter frequencies were of the order of the - 
first natural bending frequencies even though the modes were of the / 
bending-torsion type. At high angles of attack; the modes were essen- 
tially pure torsion and the frequencies were nearly equal to the first 
natural torsion frequencies. Іп most cases where stall flutter occurred, 
it occurred below maximum static lift. 


INTRODUCTION 


Flutter may be classified according to the nature of the unsteady 
flow by which it is sustained. Accordingly, that class of flutter ассош- 
panied by attached or potential flow has been termed classical fiutter. 
lt may occur in subsonic, transonic, or supersonic flow and may involve 
either single or multiple degrees of freedom. However, according to 
potential theory, classical flutter is not dependent upon angle of attack. 
А second type of flutter is associated with separated flow and is related 
to angle of attack. If flutter occurs at high angles of attack it is 
often termed stall flutter. However, a separation-type flutter may occur 
at low or moderate angles of attack at transonic speeds and it is not 
always possible to establish a definite boundary between this type of 
flutter and classical flutter. 
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The status of research on classical flutter indicates that there is 
а need for experimental research at transonic speeds. Separation-type 
flutter at high subsonic and transonic Mach numbers has become of interest 
because of its occurrence on missiles (ref. 1) and because of the need 
for data which are applicable to the design of high-speed aircraft having 
thin wings (ref. 2). The present investigation was undertaken to provide . 
such experimental flutter data for a series of rectangular wings of vary- 
ing aspect ratio. The tests were performed on the transonic bump of the 
Ames l6-foot high-speed wind tunnel at Mach numbers ranging from 0.60 to 
1.10, corresponding to а Reynolds number range from 1.1 to 2.0 million. 


SYMBOLS 
А aspect ratio, full span 
ЕТ bending stiffness, lh-in.2 
GJ torsional stiffness, lb-in.” 
Іс, polar moment of inerti& about the elastic axis per unit length, 
slug-ft?/ft 
M Mach number 
V velocity, ft/sec | 
T reduced velocity 
b wing semichord 
Е damping coefficient 
Ес gtructural Чашрїпд coefficient for the torsion mode 
£y structural damping coefficient for the bending mode 
m mass per unit length, slugs/ft 
Ig, dimensionless radius of gyration of wing section, /та/лы? 
Xc.g. center of gravity location from leading edge, percent chord 
Хаа, elastic axis location from leading edge, percent chord 
а angle of attack, deg 
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Chip minimum angle of attack at which stall flutter occurred, deg 
u mass density ratio, Z gb? 

x 
p tunnel air density, slugs/cu ft 
8, air density аб sea level, slugs/cu ft 
Woy first natural torsional frequency, radians/sec 
шы first natural bending frequency, radians/sec 

MODELS 


The seven rectangular wings tested had symmetrical NACA 64A-series 
sections with thickness-to-chord ratios of 0.02, 0.04, or 0.06. The 
models had O.5-foot chords and were constructed of solid steel or alumi- 
num. Two of the aluminum wings, 4 percent and 6 percent thick, were 
modified by cutting spanwise slots 0.005. inch wide in the upper and lower 
surfaces to reduce the torsional stiffness. The models are designated 
in this report by a number which indicates the thickness (2, 4, or 6 per- 
cent) followed by a letter to signify the material (A or 5). A letter Х 
is used to designate a slotted wing. Thus, wing ДАХ indicates а 64A00} 
aluminum wing with spanwise slots. 


The following physical characteristics are given for each wing in 
table I: mass density ratio, radius of gyration, bending stiffness, tor- 
sional stiffness, elastic axis, and center of gravity. Values of бол 
and Бос, based on the experimentally determined natural first bending and 
first torsional frequencies are given in table ЇЇ for various aspect 
ratios of each wing. The structural damping coefficients of wings 25, 
48, kA, апа АХ which were determined experimentally are presented in 
table III. 


EXPERIMENTAL DETERMINATION OF THE ELASTIC AXES AND 
NATURAL FREQUENCIES 


The elastic axis of each wing was determined by noting the movement 
of & reflected light beam while & weight near the tip was moved in the 
chordwise direction. The distance of the weight from the leading edge 
for which torsional deflection of the wing could not be detected is given 
in table I for each model in percent chord. Тһе wings were БО mounted 
as to have their meximum aspect ratios for these measurements. 
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The resonant frequencies were determined in the wind tunnel prior 
to testing. The equipment used for their determination is shown in the 
schematic and block diagrams of figure 1. The output of an audio oscil- 
lator was amplified and used to drive two 12-inch speakers. The speakers, 
mounted above and below the model, were used to excite the model through 
two strings attached to the speaker cones. The output of a crystal vibra- 
tion pickup mounted at the wing tip and a reference signal from the 
driving oscillator were used to display Lissajous patterns on an oscil- 
loscope. These patterns permitted the operator to ascertain the resonant 
frequency. 


For the purpose of exciting the bending mode, the speakers were 
located in line with the elastic axis; for exciting torsion, the speakers 
were offset to each side of the elastic axis and their polarities were so 
related that a couple was produced. 


EXPERIMENTAL DETERMINATION OF STRUCTURAL DAMPING 


The structural damping of wings 2S, 4S, 4A, and 4AX at various 
aspect ratios was determined by means of a free-oscillation technique. 
Bending oscillations were excited by the sudden release of a weight sus- 
pended from the tip at the elastic axis. А weight of 6.25 pounds мав 
used on all wings except 26 for which a weight of 5 pounds was used. 
Torsional vibrations were excited by a sharp blow struck on the leading 
edge near the tip. This procedure unavoidably introduced a small amount 
of bending which probably influenced the calculated damping coefficients. 


The motion of the wing tip was detected by a small crystal acceler- 
ometer whose mass was considered to have negligible effect on the results. 
The output of this accelerometer was recorded on an oscillograph. From 
these records the logarithmic decrement was determined for intervals of 
about 20 cycles and the structural damping coefficient was calculated 
according to the following relationship (see ref. 3): 


g = = (logarithmic decrement) 


The coefficients so determined from several records were averaged and 
the results are presented in table III. | 


TUNNEL TEST APPARATUS AND INSTRUMENTATION 


The tests were conducted on the transonic bump of the Ames 16-Ғоо% 
high-speed wind tunnel. The wings were movable spanwise through clamps 
contoured to the wing profiles. Тһе models were во positioned ав to 
provide the desired aspect ratio and were then clamped by a force of 
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800 pounds obtained through the action ОҒ an air cylinder. The aspect 
ratio was varied from 2 to 6 for the h- and 6-percent-thick wings and from 
1 to h for the 2-percent-thick wings. The maximum angle of attack was 
limited to 19.59 by the angle-of-attack mechanism. 


The flutter frequencies were determined with the aid of a modified 
phonograph crystal cartridge mounted between the upper clamping block and 
the tunnel wail. The wing vibrations were transmitted to the pickup 
through a short stiff wire resting on the upper surface of the wing. When 
the vibrations were at a constant frequency, the signal from the pickup 
was used in conjunction with a reference signal from an audio oscillator 
to display Lissajous patterns on an oscilloscope. 


А high-speed motion-picture camera (1100 frames per second) was 
mounted outside the test section, opposite and slightly above the model 
and was used to photograph the upper surface and the wing tip when flut- 
ter occurred. The motion-picture records were used to aid in the quali- 
tative determination of the flutter modes. 


METHOD OF TESTING 


The Mach number range investigated was from 0.60 to 1.10, correspond- 
ing to a Reynolds number range from 1.1 to 2.0 million. The variation of 
velocity and relative density with Mach number is given in figure 2. At 
each test Mach number selected, the model was set at O° angle of attack 
and so clamped as to provide its minimm aspect ratio. The angle of 
attack was then increased until (a) flutter was noted, (b) excessive sta- 
tie stresses or vibrations were encountered, or (с) a maximum angle of 
19.5° was reached. The model was under continuous observation in order 
to prevent its failure and the observer made notes concerning its behav- 
ior. These notes are the basis for classifying flutter conditions as 
being mild or severe. When steady flutter occurred the frequency was 
determined and in some instances a motion-picture record was taken. The 
angle of attack was then decreased to zero, the aspect ratio was increased 
by 0.33, and the procedure was repeated. Motion-picture records were 
taken of 21 cases of intermittent flutter at low angles of attack; at 
high angles of attack, 05 steady stali-flutter conditions were recorded. 


The maximum allowable angles of attack were estimated on the basis 
of static wing loads and a yield strength of 70,000 1Ъ/ва in. for steel 
and 50,000 ib/sq in. for aluminum. These limits were not rigidly adhered 
to and were subject to change during the test at the discretion of the 
observer. 


The effectiveness of the wing clamp was checked by retesting 
observed stall-flutter conditions (1.6., M, A, ap) while the wing clamp 
was rigidly bolted. The angles of attack at which stall flutter occurred 
with the wing clamp bolted were the same as those when the clamp was 
actuated by the air cylinder. 
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RESULTS AND DISCUSSION 


The results of the tests are presented in figures 3 to 9 for а11 
wings tested. Ко flutter was obtained for aspect ratios other than those 
indicated in the figures. А11 flutter data presented represent either 
Steady or intermittent flutter &s indicated by & steady or intermittent 
Lissajous pattern on the oscilloscope. 


In some instances (fig. h(&)) where mild flutter was noted at 
a = O°, the vibrations ceased with but a moderate increase in angle of 
attack and reappeared at higher angles as stall flutter. Where the flut- 
ter was more violent, particularly for the higher aspect ratios of wings 
2S and 2A (figs. 3 and 4), the vibrations persisted to higher angles of 
attack and limited the investigation to the angles shown. It was not . 
possible to establish a definite boundary between classical and separation- 
type flutter in all instances. Only the steady flutter (figs. 3 to 9) at 
high angles of attack and the flutter which occurred initially at low 
angles (09 to 29) were considered to be subject to classification as 
stall flutter and low-angle-of-attack flutter, respectively. 


Flutter at Low Angles of Attack 


The determination of specific flutter points at low angles was not 
as decisive as was the case at high angles. The motion pictures taken 
during the tests indicated that the flutter of wings 2А and 2S at small 
angles of attack was of the coupled bending-torsion type in some instances; 
however, the predominant frequency was of the order of the first natural 
bending frequency. The vibrations of wing 2A (fig. 4) were quite severe 
from Mach number 0.80 to 1.10 for aspect ratios of 2.33 and 2.67. Тһе 
vibrations were very violent and of large amplitude for aspect ratio 3.00 
and the predominant frequency from Mach number 0.85 to 1.00 was approxi- 
mately 50 cycles per second, The model failed at M = 1.02, A = 3.00 
ав the tunnel speed was being reduced. The failure is believed to have 
been а result of structural fatigue rather than of divergent flutter since 
the model appeared to be more stable at Mach numbers above 1.00. 


The intermittent flutter noted for wing 2S (fig. 3) was near the 
first natural bending frequency for aspect ratios 2.67 to 3.33. Те 
vibrations for aspect ratios 3.67 and 4.00 became progressively more 
severe with increasing Mach number, so that the Mach number was limited 
to 0.94. No flutter frequencies were determined for these aspect ratios. 


The vibrations of wing 4AX (fig. 7) at low angles were also near the 
first natural bending frequency for aspect ratios of 3.33, 3.67, and h.00 
and did not appear to be in a coupled mode at any time. The vibrations 
of this wing when positioned at aspect ratios greater than 4.00 were too 
erratic to determine the frequency. 
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No low-angle-of-attack flutter was encountered for wings 4S, hA, 
6A, or бАХ. Тһе intermittent flutter at moderate angles of attack which 
was noted for wing 4A (fig. 6; A = 3.67, 4.00; а = 3.59) and wing 6A 
(fig. 8; A = 6.00; а = 4.5°) was not considered to be of the low-angle- 
of-attack type since no flutter occurred at a = 09, 


Flutter at High Angles of Attack 


The stall flutter was of а steady torsional mode at or near the 
first natural torsional frequency of the wing, a phenomenon which sub- 
stantiates results of previous investigations. The predominant torsional 
mode was clearly evident in the motion pictures for various stall-flutter 
conditions. At a given Mach number, the intensity of the vibrations 
generally increased with angle of attack to the limits of the test. How- 
ever, for wing 25 at А = 1.67 (fig. 3), it was possible to increase the 
angle of attack until a stable condition was again established. 


Several of the wings tested appeared ta be less susceptible to stall 
flutter at the higher test Mach numbers. This trend was particularly 
evident in the case of wing 4AX (fig. 7; A = 2.00) for which steady flut- 
ter occurred at M = 0.90 (о = 16.57), latem en flutter occurred at 

= 0.94 (a = 15.59 to 19.59), and no flutter was noted from M = 0.98 
to М = 1.10 at angles of attack up to 19.59. A similar but less pro- 
nounced trend was noted for this wing at А = 2.33 апа 2.67, as well ав 
for wing бАХ (fig. 9; А = 3.33 to 4.67), and wing 28 (fig. 3; A = 1.67). 


Ав mentioned previously, the occurrence of steady stall flutter was 
pronounced and easily determined by the Lissajous patterns on the oscil- 
loseope. The minimum angle of attack at which steady flutter occurred 
(ағ) has, therefore, been plotted &s & function of Mach number in figures 
lota) to 15(а) to summarize the effects of Mach number and aspect ratio 
for all wings except 6A. This wing exhibited only a mild form of inter- 
mittent flutter within the angle-of-attack range investigated. 


Data from part (a) of figures 10 to 15 are presented in & different 
form in part (b) of these figures to introduce the conventional flutter 
parameter V /büq, ав а function of ар. Тһе parameter Ч/рес, decreased 
rapidly with increasing angle of attack. А value of V/bug 2 1 has been 
suggested in reference 2 &8 being useful in estimating the minimum stall- 
flutter velocity of thin wings. Explicit verification of this criterion 
is not possible, however, since no date were taken at Mach numbers less 
than 0.60 or at angles greater than 19.59. 


In figure 16, ap is plotted as a function of aspect ratio at 
M = 0.60 to indicate the effect of а reduction іп torsional stiffness 
on оғ for the h-percent and Ó6-percent aluminum wings. Ав noted previ- 
ously, wing 6A exhibited only mild intermittent flutter up to the maximum 
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angles of attack allowable from mechanical or strength considerations. 
These values of о, shown in figure 16 for this wing, аге considered to 
be conservative estimates of the true values of ар. The data ss pre- 
gented indicate that reduction of torsional stiffness had a more pro- 
nounced effect for the 6-percent wing than for the 4-percent wing. А 
similar comparison at other Mach numbers is not possible, however, because 
of the lack of data for the solid wings. Wing 4A failed during the test 
before sufficient siall-flutter data were obtained. The data for wing 

6A at Mach numbers higher than 0.60 were not deemed sufficiently reliable 
to continue the comparison. 


Stall flutter has been associated in the past with the static lift 
curve (see ref. 4). This association led the authors of reference 4 to 
suggest three possible causes of stall flutter: (1) static instability 
due to the negative slope of the lift curve &t angles beyond the stall, 
(2) dynamic instability due to hysteresis at the stall, end (3) insta- 
bility resulting from Kármán vortex excitation. 


In figures 17, 18, and 19 the angles of attack at which flutter was 
Observed have been noted on static lift curves taken from reference 5. 
These static data were obtained from tests conducted in the same facility 
апа on models of the ваше size ав those of the subject investigation. Тї 
must be noted, however, that these curves are for wings having symmetrical 
NACA 63A-series sections rather then the 64A-series used in the present 
investigation. Comparison of availeble 63-series and 6h-series-airfoil- 
section data in references 6, 7, and 8 indicates that the static data 
used are indicative of the true static characteristics. 


The limited number of available data ав shown in figures 17 to 19 
indicate that with respect to the three possible causes of stall flutter 
cited above, only the steady flutter of wing 4AK at the lower aspect ratios 
(fig. 18) and the intermittent flutter of wing 6A, A = 6.00 (fig. 19) may 
be attributed to static instability due to the negative slope of the lift 
curve. The remaining cases of steady flutter can then be considered as 
probably being a result of the hysteresis loop described by the lift curve 
since the stall-flutter frequencies were considerably lower than the vor- 
tex frequencies calculated by the method given in reference 4. 


CONCLUDING REMARKS 


Flutter at low angles of attack occurred only for the 2-percent- 
thick steel and aluminum wings and for а 4-percent-thick aluminum wing 
slotted spanwise to reduce the torsional stiffness. The predominant flut- 
ter frequency was of the order of the first natural bending frequency. 


Stall flutter occurred at frequencies which were in good agreement 
with the experimentally determined first natural torsional frequencies. 
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The angles of attack at which stall flutter occurred have been 


indicated on the static lift curves of а series of rectangular wings 
having symmetricai НАСА 63A-series sections. Іп the majority of cases, 
flutter occurred ал angles below those for maximum lift. 


Ames Aeronautical Laboratory 


National Advisory Committee for Aeronautics 
Moffett Field, Calif., Jan. 29, 1954 
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TABLE I.- CHARACTERISTICS ОР THE WINGS TESTED 


OT 


Torsional | Хе.в., | Хс.р.; 
BtiffnesB,| stiffness, 


ЕТ, GJ, 
lb-in.? 1Һ-іп.2 


1Spanwise slots having a depth of 55 percent of the wing-section ordinate located 5 
at chordwise Intervals of 0.3 inch between 5- anā 90-percent chord. 

“Spanwise slots having a depth of 35 percent of the wing-section ordinate located 
at intervals of 0.3 inch between 5- and 90-percent chord. 

Note: Values of mass density ratio ы are based on the density of air at sea 
level. For variation of test density with Mach number see figure 2. 
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TABLE II.- STILL-AIR FORCED-VIBRATION CHARACTERISTICS OF WINGS TESTED 


+ ONA 


нн 
(ШТ 
Со = CO ON TO afr 


lb wing semichord, 0.25 fi. 
Suh first natural bending frequency, radians/sec. 
Sw, first natural torsional frequency, radians/sec. 


TABLE III.- EXPERIMENTALLY DETERMINED STRUCTURAL DAMPING COEFFICIENTS 


Aspect ratio 


16, structural damping coefficient for the bending mode. SMA 
ка structural demping coefficient for the torsion mode. 
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Figure 1.- Sehematic and block diagrams of apparatus used for the 
determination of the natural frequencies. 
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Figure 2.- Variation of velocity and relative density with Mach number. 
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Figure 3.- Flutter observations, wing 25. 
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Flutter observations, wing 2А, 
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Figure 4.- Concluded. 
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Figure 5.- Flutter observations, wing 49. 
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Figure 5.- Concluded, 
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Figure 6.~ Flutter observations, wing 4A. 
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Figure 6.- Concluded. 
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Figure Т.- Flutter observations, wing HAX. 
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(а) Variation of ар with Mach number. 
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Figure 10.- Effect of aspect ratio and Mach number on the stall flutter 
of wing 25. 
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Figure 11.- Effect of aspect ratio and Mach number on the stall flutter 
“ОҒ wing 2А. 
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Figure 14.- Effect of aspect ratio and Mach number on the stall flutter 
of wing БАХ. 
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Figure 15.- Effect of aspect ratio and Mach number on the stall flutter 
of wing бАХ. 
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Figure 16.- Effect of a reduction in torsional stiffness en stall flutter 
at Mach number 0.60. 
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Figure 17.- Stall flutter of the 644002 wings related to the static lift 
curves of а 63A002 wing. 
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Figure 18.- Stall flutter of the GhAO0h wings related to the static lift curves of а 63A004 wing. 
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lift curves of a 63A006 wing. 


~ Stall flutter of the 614006 wings related to the static 
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